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I— I We present results for the propagators and the ghost-gluon vertex of Landau gauge Yang-Mills 

"5 theory obtained from Dyson-Schwinger equations. Solving these three quantities simultaneously 

pi I constitutes a new step in truncating these equations. We also introduce a new model for the three- 

gluon vertex that is motivated by lattice results. It features a zero crossing which is confirmed 
a posteriori by a Dyson-Schwinger calculation. Within our setup we can reproduce lattice data 
,_l very well. We establish that also for the ghost-gluon vertex a difference between decoupling 

and scaling solutions is present. For the scaling solution we discuss the possibility of modifying 
QQ the infrared exponents via an angle dependence of the ghost-gluon vertex. However, no such 

dependence is found in our calculations. Finally, we calculate the Schwinger function for the 
gluon propagator. 
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1. Two-point functions 

Extracting non-perturbative information from a quantum field theory is in most cases a hard 
task. One possibility is to calculate its Green functions which are related by several sets of func- 
tional equations. A concrete example of such a set are the equations of motion of the Green func- 
tions, called Dyson-Schwinger equations (DSEs) [1, 2]. An inherent problem is that there ai^e 
infinitely many equations. The usual way to make quantitative calculations is to neglect some 
Green functions and/or model others. The remaining ones are then solved self-consistently. The 
natural way to test the reliability of such a truncated system is to compare it with a larger trunca- 
tion. However, in practice this is a difficult task and it can take some time until such a calculations 
becomes feasible both from the conceptual and from the computational point of view. An example 
for this is provided below. Nevertheless, it has to be stressed that raising the number of included 
Green functions can allow a systematic improvement. In practice, however, it is more common to 
estimate truncation artifacts from comparisons with other methods. 

Sometimes we are in an even better situation and can estimate the influence of truncations. 
The most prominent example is surely the high energy regime of QCD. Since the strong interaction 
is asymptotically free, a perturbative loop expansion provides a reliable approximation which can 
be improved systematically by including higher powers of the coupling. On the other end of the 
energy spectrum, in the infrared (IR), it can also be possible to determine a clear hierarchy of 
leading and suppressed contributions, if a so-called scaling solution is realized [3, 4, 5]. Then the 
qualitative behavior of all Green functions can be determined in form of IR power laws. Landau 
gauge Yang-Mills theory may actually possess such a solution [6, 7, 8, 9]. Knowledge of additional 
information like this is of course a useful guideline in constructing truncations. 

In the present work we consider Landau gauge Yang-Mills theory. It is illustrative in the light 
of the discussion above to recall the development of truncations in this case. In a very basic approx- 
imation the ghosts and all four-point functions were neglected [10]. The only DSE to be solved 
is the gluon propagator DSE. The resulting gluon propagator is IR divergent, which is nowadays 
known to be wrong. The next step is the inclusion of ghosts, so that a coupled system of two 
DSEs has to be solved. This step was taken 18 years later [6, 7] and resulted in a drastic change, 
as the gluon dressing became IR suppressed. Based on analytic and numeric results from various 
methods, we can be confident that a truncation at this level yields quaUtatively reliable results. 

Within this class of truncations two types of solutions can be obtained. However, this is not 
related to the truncation itself but to the boundary conditions imposed on the DSEs and a difference 
is only visible at small momenta [8]. The decoupling type consists of a family of solutions which 
have an IR finite ghost dressing and an IR finite gluon propagator [8, 11, 12, 13, 14]. The scaling 
solution, on the other hand, features an IR divergent ghost dressing and an IR vanishing gluon 
propagator [6, 7, 8, 9]. It is still under debate if both solutions exist, especially since all "standard" 
lattice calculations recover the decoupling type of solutions. However, lattice results do depend 
on the method used for sampling Gribov copies, see, e.g., [15, 16]. Recently a way was found to 
sample Gribov copies such that the results are modified as expected from functional equations [17]. 

In the following we will present results from a new truncation that goes beyond the propagator 
level by including the ghost gluon-vertex dynamically and employs a model for the three-gluon 
vertex in agreement with lattice results [9]. 
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2. Three-gluon vertex 

For the three-gluon vertex we employ the following model motivated by lattice calculations: 

(2.1) 

with the damping factors 

:= {1.1) 

We choose hjR <0 so that the vertex becomes negative for low momenta. This is observed with 
Monte Carlo simulations in two and three dimensions [18, 19], but not in four [19], probably 
because the lattice data does not reach far enough into the IR. We performed a DSE calculation for 
the vertex based on the results obtained for the propagators and the ghost-gluon vertex and taking 
into account only the IR leading diagram. We indeed found a zero crossing, see Fig. 1, which is 
compatible with the available lattice data. The position of the crossing may still change in a more 
advanced truncation, but it is very unlikely that it will vanish. This is con^oborated by results in two 
dimensions, where more detailed calculations were performed and the influence of other diagrams 
was found to be negligible in the IR [20]. 

In Eq. (2.1) the IR leading part is damped out for larger momenta by the damping functions 
f{x). The other term is motivated by perturbation theory: By choosing a and j8 appropriately, we 
reproduce the anomalous dimension of the three-gluon vertex, which is 73^ = 1 + 35 = 17/44 with 
5 = —9/44 the ghost anomalous dimension. The second condition to fix a and j3 is IR finiteness. 
We obtain for the scaling solution a = —1 — 65 and j8 = — 1 — 35 and for decoupling a = 3 + 1/5 
and j3 = 0. In order to produce the correct anomalous dimensions for the propagators, it is necessary 
to make an adjustment to the gluon loop in the gluon DSE, where the renormalization constant Zi 
of the three-gluon vertex appears. We replace it by a momentum-dependent renormalization group 
improvement term that equals the UV part of the three-gluon vertex model [7]. Such a modification 
is unavoidable in order to obtain the correct anomalous dimensions, because they result from the 
resummation of infinitely many diagrams which cannot be taken into account directly in the present 
approach. 

In Fig. 1 results for the gluon dressing function for varying values of A^g are shown. The 
influence of the zero crossing is clearly visible: If its position moves to higher momenta, the bump 
in the gluon dressing function rises. In general the model shows good agreement with the lattice 
data, which is, however, for two colors while our calculation is for three colors. From Fig. 1 we 
conclude that the zero crossing, if chosen at a value in agreement with lattice data, has no influence 
on the propagators. The reason is that in the gluon loop the vertex is multiplied with two gluon 
propagators which are IR suppressed. Besides the IR part also the Bose symmetry of the vertex 
has an effect on the mid-momentum regime of the propagators, see ref. [9]. In the following we 
will use the three-gluon vertex model to mimic effects of the two-loop diagrams in the gluon DSE 
by choosing the zero crossing such that the propagators agree as good as possible with lattice data. 
The model with the corresponding choice of parameters is called optimized effective vertex. 
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Figure 1: Comparison of different parameters for the three-gluon vertex model: hjR = — 1, Aa^ = 
0, 0.86, 1.08, 1.31, 1.53, 1.1 GeV . Left: Symmetric configuration of the vertex model. Thick lines from left 
to right correspond to rising values of the parameter Asg. The thin line is the IR leading contribution of the 
three-gluon vertex calculated from its DSE showing a zero crossing at very low momenta. At zero momen- 
tum it is finite. Individual points are lattice data from ref. [19]; Red circles correspond to N = l6/j5 = 2.2, 
green squares toN = 16//3 = 2.5, blue diamonds toN ^ ll/ji — 2.2 and orange triangles to N = 22/)3 =2.5. 
Right: The corresponding gluon dressing function. Curves from bottom to top correspond to values of A3g 
from low to high. 

3. Ghost-gluon vertex 

The ghost-gluon vertex is an important quantity in modern truncation schemes of functional 
equations, because it is the only primitively divergent vertex via which the ghost interacts. The 
reason for the success of modem truncation schemes is that this vertex has only a very mild mo- 
mentum dependence. Thus approximating it by its bare expression produces reasonable results. 
Non-perturbative modifications of the vertex were investigated with various approaches: Monte- 
Carlo simulations [15, 18, 19], semi-perturbative DSE calculations [21, 22], OPE analysis [23] and 
based on this indirect studies via effects in the ghost DSE [24]. Here we will include the vertex 
dynamically, i.e., we solve the DSEs of the two propagators and the vertex simultaneously. We 
employ the following basis for the vertex: 

lf''^'"''{t,p,q) := f"Tf,''{k-p,q) := igf'^' {p^A{t,p,q)+k,B{k-p,q)) , (3.1) 

where the momenta k, p and q refer to the gluon, the anti-ghost and the ghost, respectively. Due 
to the transverse gluon propagator of Landau gauge, only the dressing A{k;p,q) is relevant. The 
corresponding DSE as well as those for the propagators can be found in ref. [9]. 

The effect of the dynamic ghost-gluon vertex is illustrated in Fig. 2, where calculations with a 
bare ghost-gluon vertex are compared to one with a dynamic vertex. The red/continuous line cor- 
responds to the solution with the optimized effective three-gluon vertex and a dynamic ghost-gluon 
vertex. The coiTcsponding calculation with the same parameters for the three-gluon vertex but a 
bare ghost-gluon vertex is represented by the blue/dotted line. However, although there is quite 
some change in the mid-momentum regime of the gluon propagator, this only provides a quali- 
tative picture, since the model for the three-gluon vertex depends on the propagators and is thus 
different in the two calculations. Finally, we note that it is possible to approximately accommodate 
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Figure 2: Comparison of the gluon propagator dressing function Z{p^) and the ghost dressing function 
G{p^) for different setups. Red/continuous Une: dynamic ghost-gluon vertex, = 2.1 GeV. Green/dashed 
hne: bare ghost-gluon vertex, Aa^ — 1.8 GeV. Blue/dotted line: bare ghost-gluon vertex, A^g = 2.1 GeV. 



the modifications introduced by a dynamic vertex also in tlie tliree-gluon vertex model. Tliis can 
be seen from the green/dashed line, for which the parameters of the model were adjusted such that 
the gluon propagator from the lattice is reproduced using a bare ghost-gluon vertex. 

4. Results 

Technical details on the calculations can be found in ref. [9]. The programs DoFun [2, 25] and 
CrasyDSE [26] were used for deriving and solving the DSEs. 

The resulting dressing function for the propagators are shown in Fig. 3. As can be seen, it is 
possible to choose the parameter A^g such that we obtain a very good agreement with lattice data. 
We would like to stress that this is true for both dressings. The Schwinger function A{t) of the 
gluon propagator is shown in Fig. 4. It turns negative at about 2.2 fm. This value is higher than 
observed in previous calculations [8]. Furthermore, a second zero crossing is observed. Since the 
Schwinger function is partially negative, positivity is violated and the gluon cannot represent an 
asymptotic state of the physical state space. 

Results for the ghost-gluon vertex are shown in Fig. 5. The depicted momentum configura- 
tion corresponds to almost anti-parallel gluon and anti-ghost momenta. As expected the dressing 
is constant in the IR and the UV and has a small bump in the intermediate momentum regime. 
The ambiguity between scaling and decoupling is also present in the vertex, as the IR constant is 
different: For decoupling the dressing is 1, because the loop corrections are IR suppressed by the 
gluon propagators. For scaling, however, they yield a non-vanishing contribution which lifts the 
dressing above 1. This can also be seen from the red/continuous line in Fig. 4, where additionally 
two different decoupling solutions are plotted. 

For the scaling solution the exponents of the IR power laws for all dressing functions can be 
calculated in terms of one variable K [3, 4], which can be calculated analytically [29, 30]. Under 
the assumption that the IR limit of the ghost-gluon vertex is regular, its value is k: = 0.595353. This 
value can only be modified by a non-analyticity at zero momentum [30]. One possibility for this, 
a dependence on the angle between two momenta, was studied in ref. [20] . Here we tried to find 
a deviation by not fixing the value of K manually, as is normally done, but leaving it free. Thus 
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Figure 3: Comparison of the results for the gluon propagator dressing function Z{p^) and the ghost dressing 
function G{p^) with lattice data [15]. The red/continuous line was obtained with a dynamic ghost-gluon 
vertex and the optimized effective three-gluon vertex with hm = — 1 and As,, =2.1 GeV . For reference the 
green/dashed line is shown, which was obtained with the vertex models of refs. [27, 28] and a bare ghost- 
gluon vertex. Lattice data [ 1 5] is for j3 = 6 and lattice sizes of L = 32 (blue circles) and L = 48 (orange 
squares). 

|A(t)| 




Figure 4: Left: The absolute value of the Schwinger function of the gluon propagator. In the region between 
approximately 2.2 GeV and 9 GeV the function is negative. Right: The ghost-gluon vertex at the symmetric 
point. The red/continuous line is from the scaling solution, and the green/dashed and blue/dotted lines from 
two different decoupling solutions. 

the ghost-gluon vertex could develop an angle dependence, which was, however, not observed. 
The differences between the vertices as obtained from both procedures were below 0.1% and thus 
within numerical accuracy. The value of K was found to be a little bit higher, slightly above 0.6. 
Thus we did not find any non-analyticities in the ghost-gluon vertex. 

5. Summary and conclusions 

In general DSEs are a useful tool for non-perturbative investigations. For Landau gauge Yang- 
Mills theory we shortly recalled the history of different truncations and explicitly demonstrated that 
further improving the truncation leads only to quantitative effects. One improvement presented here 
is the dynamical inclusion of the ghost-gluon vertex. We find that the ambiguity of decoupling and 



6 



Going beyond the propagators of Landau gauge Yang-Mills theory 



Markus Q. Huber 




scaling solutions persists also at the level of three-point functions. We also introduced a model 
for the three-gluon vertex that properly describes lattice data. This model may also be relevant for 
future calculations when the so-called squint diagram is included into the calculation. Results for 
the sunset were presented in [31]. In summary, we could reproduce both propagators as obtained 
from Monte Caiio simulations within the present setup. 
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